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Abstract. A homogeneous sample of galaxy redshifts in 
the core regions (R < 0.5 h" 1 Mpc) of 12 clusters is used to 
measure the frequency of substructure with different tests. 
In 50 % of the cases substructure is detected, a frequency 
which agrees well with previous studies of cluster cores. 
Magnitude information and rough morphological classes 
are also available for 80 % of the sample, allowing us to 
confirm that bright galaxies (M< -22 ma9 ) have a signif- 
icantly lower velocity dispersion than the rest. Elliptical 
galaxies are the main responsibles for this luminosity seg- 
regation in velocity space, whereas no segregation can be 
seen for spiral galaxies. Given the coincidence of substruc- 
ture and luminosity segregation, a cluster model with an 
old population of ellipticals which are under the effect of 
dynamical friction in each subcluster is thus favoured by 
these observations. Spiral galaxies seem to be late arrivals, 
or are passing in front or behind the core of the cluster. 

Key words: Galaxies: clusters: general 



1. Introduction 

One of the main results of recent galaxy cluster research 
is the unambiguous finding that these high-density struc- 
tures are still continuing to be built or at least significantly 
reshaped. Undoubtedly, the discovery of substructures has 
played an important role, changing our understanding 
about the degree of dynamical evolution in galaxy clus- 
ters. 

Several statistical treatments to measure the frequency of 
substructures have been presented during the last decade. 
Nevertheless, the question about the significance of sub- 
structure detections on small scales or in the core regions 
of clusters has been touched only in a few works (Fitchett 
& Webster 1987; Mellier etal. 1988; Escalera etal. 1992; 
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Salvador-Sole etal. 1993). 

It is interesting that the frequency and degree of dumpi- 
ness in the centers of clusters could be helpful in estab- 
lishing the density profile of dark matter, and even allow 
an estimation of fl (e.g. Richstone etal. 1992). If namely 
dark matter was strongly concentrated towards the clus- 
ter center one would expect tidal forces acting towards the 
disruption of subclumps (see Gonzalez-Casado et al. 1994 
for a discussion). 

Examining the very centers of galaxy clusters is rewarding 
for another reason. The importance of accretion or "canni- 
balism" (Ostriker & Hausman 1977) during the formation 
of a cluster is not well established, because right in the 
cores of galaxy clusters the galaxy velocity dispersion is 
expected to be too high to allow an efficient dynamical 
friction (Merrifield & Kent 1989, Gebhardt & Beers 1991, 
Blakeslee & Tonry 1992). An attractive solution to this 
problem has been proposed by Merritt (1985). 
As also numerical simulations suggest (West & Richstone 
1988, Serna etal. 1994), dynamical friction could be an 
important mechanism during the evolution of galaxy clus- 
ters. If this is the case, it should be possible to detect 
signs of equipartition of kinetic energies for the most mas- 
sive galaxies in dense regions, leading to a correlation of 
velocity dispersion with galaxy luminosity. Early observa- 
tional studies of mass segregation were mostly limited to 
the positions of galaxies in projection, with the exception 
of e.g. Chincarini & Rood (1977), and no general agree- 
ment was foud about the relevance of luminosity segrega- 
tion in velocity space. Recently, Biviano etal. (1992) once 
again stated the evidence for mass segregation in the ve- 
locity distribution of a merged cluster sample taken from 
the literature. Presumably, the phenomenon is easily over- 
looked in individual clusters studies because of the limited 
number of galaxies involved. 

On the other side, a difference in the kinematic properties 
between early and late-type galaxies has also been de- 
tected (Binggeli etal. 1987, Sodre etal. 1989). This means 
that one should look separately at early and late-type 
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galaxies, because galaxies which didn't take part in most 
of the cluster evolution are not expected to show signs of 
mass segregation. 

This paper addresses the problem of mass segregation and 
central substructure based on a homogeneous sample of 
kincmatical data (Stein 1996, hereafter Paper I). It is or- 
ganized as follows: Section 2 gives a description of the data 
used, i.e. the source of redshifts, magnitudes and types for 
the galaxies, as well as cluster definition and selection. The 
analysis of substructures on individual clusters is done in 
Section 3, while luminosity segregation is measured on the 
merged cluster sample (Section 4). The results are then 
discussed in Section 5. 

2. The sample 

The present work is based on a rcdshift catalogue involv- 
ing the central regions of 15 clusters of different richnesses 
in the redshift range 0.01 < z < 0.06 (Paper I). Note that 
several of the clusters are well studied, nearby objects 
like Centaurus and Abell 400 or have been extensively 
observed for specific purposes (e.g. Abell 3558 : Bardclli 
etal. 1994). Inside these regions nearly complete (80 %) 
rcdshift information down to Ms w 18 ma9 could be ob- 
tained, as well as raw morphological types for most of the 
galaxies. The present sample differs from those of previ- 
ous investigations because it concentrates onto the very 
cores of the clusters (R < 0.5 h _1 Mpc), and because of 
the homogeneity of kinematical information. 
Cluster centers were chosen in order of preference from (a) 
published X-ray centers (Lahav etal. 1989, Edge & Stew- 
art 1991), (b) the position of the cd-galaxy, if present, (c) 
the Abell cluster center or (d) by taking the center of the 
relevant Optopus field (see Paper I), which roughly coin- 
cides with the galaxy density peak. The only exception to 
this procedure was Abell 1736, whose OPTOPUS field was 
offset by 0.05 h" 1 Mpc with respect to the X-ray center. 

A common selection radius of 0.5 h _1 Mpc around 
each cluster centre was therefore chosen, using approxi- 
mate mean cluster velocities and assuming all clusters to 
be located in the Hubble flow. For Hydra and Centaurus 
(z w 0.01) two adjacent fields had been taken, and for the 
other very nearby clusters data from the literature was 
available to help fill the region. 

Table 1 lists the fields and also gives the name of each 
cluster around whose centre the selection has been done 
(2), kind of center determination (3), where X : X-rays, 
cd : cd-galaxy position, A : Abell center, O : center of Op- 
topus field, as well as centre coordinates (4), and selection 
radius in arcminutes (5). 

2.1. Membership selection 

As a first step towards meaningful cluster definitions, we 
then looked at velocity histograms in the direction of all 
15 fields, eliminating obvious (5 a) fore- or background 



galaxies. In one case (field 12) two distinct structures could 
be identified : the main cluster at v = 4500 kms -1 and a 
smaller group at v = 11000 kms . Moreover, some non- 
gaussian velocity histograms with indication of bimodal- 
ity could be recognized. A statistical test was thus ap- 
plied which returns the likelihood that the biggest gap in 
the dataset could occur in a normal distribution (adapted 
from the ROSTAT package, Beers etal. 1990). Only in the 
case of field 9 a gap was found whose size was inconsistent 
at the 5 % level with an underlying gaussian distribu- 
tion. Field 9 was thus considered bimodal and subdivided 
at v = 11500 kms -1 , the big gap location. To assess fi- 
nal cluster membership, a 3-ct clipping technique (Yahil 
& Vidal 1977) was then applied to each of the 17 galaxy 
units recognized so far. Table 2 lists the number of mem- 
ber galaxies after 3-cr clipping (column 3), as well as the 
resulting limits in redshift space (4). 

Mean cluster velocities and velocity dispersions were 
then determined for each cluster in the sample, using bi- 
weight estimators (Beers etal. 1990). These estimators 
have the advantage of being robust against outliers and 
are particularly useful when working with small datasets. 
Cosmological effects are taken into consideration follow- 
ing Danese etal. (1980). The resulting values and their 
uncertainties are listed in Table 2, columns 5 and 6. 

Our cluster mean velocities and velocity dispersions 
do agree quite well with those quoted by Struble & Rood 
(1991) and Girardi etal. (1993). Nevertheless, there are 
some cases where a discrepancy is found. The reason here- 
fore must be related to the effect of substructures in the 
central cluster regions, which tend to inflate velocity dis- 
persions, and to the presence of luminous galaxies of low 
dispersion in the cluster cores (as shown in the next sec- 
tions). Nevertheless, given that our results are virtually 
free from larger scale contamination, we consider them 
good estimators of central mean velocity and velocity dis- 
persion. With the present data it was not possible to dis- 
cern the bimodal structure of Centaurus using a gapper 
test. For this reason we give a global value for the mean 
cluster velocity and for the velocity dispersion. Neverthe- 
less, substructure is detected in this and in some of the 
other clusters, as we will see below. It may be, therefore, 
that the kinematic parameters given in Table 2 do not re- 
flect the true dynamical state in 50 % of the cases. They 
should be treated as first estimations and were calculated 
because a measure for the cluster dispersion is needed for 
the Lee-test simulations (see below). 

Some of the galaxy samples are clearly too poor to 
be used for the substructure analysis and have to be dis- 
carded. Computation of the Dressier & Shectman (1988) 
substructure test requires the evaluation of velocity dis- 
persions from samples of >/N neighbours around each 
galaxy, where N is the total number of galaxies (Bird 
1994). Given that 5 is a lower limit for the determination 
of standard deviations, we chose a number of TV > 5 2 = 25 
galaxies as the minimum richness. 
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Table 1. Cluster sample 
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We are thus left with 12 clusters (flag "Y" in Table 
2, column 7) and a total of 576 galaxy redshifts, 2/3 of 
which are taken from Paper I and have mean errors well 
below 50 km s _1 . 



2.2. Photometric data 

As we wanted to look at dynamical friction effects, lu- 
minosity information was necessary. For three of the 
clusters detailed photometrical studies could already be 
found in the literature : A 3526 = Centaurus by Dick- 
ens etal. (1986), A S 805 by Millington & Peach (1989) 
and A 4038 by Green etal. (1990). In addition, bj magni- 
tudes from the COSMOS catalogue were kindly provided 
by H. MacGillivray (1993) for ten clusters of high galactic 
latitude. Given the selection of high latitude clusters no 
correction for reddening or extinction was applied. 
A comparison of COSMOS magnitudes with those of 
Green etal. (1990) in A 4038 shows an excellent agree- 
ment, with differences of only a few hundredths magni- 
tudes. For some of the brightest spiral galaxies COSMOS 
magnitudes were not available (due to the problem of re- 
solved HII regions), which meant that magnitudes had 
to be taken from a bright galaxy catalogue. It should be 
noted that magnitudes from the COSMOS catalogue are 
bj magnitudes, while those for Centaurus are G26.5 mag- 
nitudes, and for most of the other clusters generic optical 
magnitudes were taken from different sources. Magnitudes 
were not transformed to a common scale, because for our 
statistical analysis a spread of a few tenth of magnitudes 
for individual galaxies could be taken into account. 



2.3. Morphological data 

Galaxy types were searched for in the literature, mainly 
resulting in a sample from Dressler's (1980) and UGC (Nil- 
son 1973) catalogues, as well as from Huchra's (1991) col- 
lection, as implemented in the DIRA2 database (Astronet 
Data Base Group Italy, Bologna). All galaxies were then 
divided into three classes: E, SO and S. In addition, many 
of the remaining galaxies were classified by the author 
into one of the above classes, after visual inspection of 
ESO-Schmidt plates. Because careful classification is dif- 
ficult on this kind of plates, a check was made involving 
105 galaxies with independently known types. Of these, 15 
were classified by the author with an "earlier" type and 13 
with a "later" type than literature values, corresponding 
to an agreement of around 75 %. 

3. Substructures 

Ideally, substructure tests should be applied to galaxy 
samples which are complete in magnitude. Here, we will 
use all the available data without regard to the magni- 
tudes, for different reasons. First, the small sample size 
makes any further restriction unreasonable. In addition, 
complete magnitude information is not available for all 
cluster fields. Nevertheless, in those fields where magni- 
tudes were available from COSMOS, rcdshift completeness 
amounts to ca. 75 % of the galaxies down to 18 maff . Given 
the small field of view and the extensive redshift cover- 
age, we will concentrate on substructure analysis methods 
which make use of the velocity information. We chose the 
test of Dressier & Shcctman (1988, hereafter DS-test), the 
Lee test (Fitchett 1988), and several tests that check for 
departures from normality of the galaxy velocity distribu- 
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Table 2. Redshift ranges and kinematical parameters 
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tion. Note that: 

a) For the DS-test we followed the same procedure as Bird 
(1994), using VN galaxies to define the neighbourhood of 
each galaxy, where N is the total number of galaxies in- 
volved in the analysis. 

b) To assess the significance of substructure detections for 
the Lee test we performed the same test on 100 cluster 
simulations. Modelling of the simulated clusters was done 
with an 1/r galaxy density profile (Fitchett & Webster 
1987) and with the same velocity dispersion as in Table 
2. Steps of 4 degrees have been used for the orientation 
of the projection axis in the X-Y-cz plane. Note that the 
simulations for field 9b were done with an offset of 0.05 
h -1 Mpc between the peak of the 1/r density distribution 
and the center of the circular selection region. A higher 
likelihood for the presence of substructure would result if 
the offset was not accounted for (92 % instead of 87 %). 

c) Three different tests for normality have been employed, 
i.e. the U 2 , W 2 and A 2 tests taken from the ROSTAT 
package (Beers etal. 1990), the resulting significance be- 
ing the mean from these three tests, which agree very well. 
Each of these substructure tests has its strengths and 
weaknesses. For example, the DS-test measures the de- 
viation of the local from the global kinematics, allowing 
to detect efficiently small offcenter groups, but sometimes 
failing in cases where two clumps of equal size and different 
mean velocities overlap in projection. On the other side, 
the Lee test is designed to be a very general maximum- 
likelihood method, but computational constraints confine 
its use to detecting bimodal structures of comparable size, 



thus being unable to state about the presence or absence 
of multimodal structures. It is also trivial that testing the 
gaussianity of a velocity sample alone cannot give clues 
about substructures which have same means and disper- 
sions but differing locations in the plane of the sky. We see 
that all of these tests are bound to miss some manifesta- 
tions of substructure, each one being sensitive to some par- 
ticular configuration. For these reasons, a combination of 
all methods should allow a better judgement to be made. 
Our statement about the existence of substructure relies 
upon the fact that at least one of the three detection meth- 
ods could find significant signs for it. Results are shown 
in Table 3, which gives field number (column 1), name of 
the cluster (2), significance level for the tests mentioned 
above (3-5), and total significance level for substructure 
(6), which is the highest value of columns 3, 4, and 5. As 
can be seen from the analysis results, 50 % of the clusters 
show clear signs of substructure (5 % significance level). 
This frequency of substructure in cluster cores has to be 
considered a lower level because of the intrinsic property of 
individual tests to miss some manifestations of substruc- 
ture. It should be noted that the same order of magnitude 
for substructure frequency has been found by other au- 
thors (e.g. Escalera etal. 1994; West 1994), in particular 
also involving comparable spatial resolution but without 
redshift information (Salvador-Sole etal. 1993). 

4. Segregation in velocity 

Radial velocities and magnitudes had to be normalized 
before a kinematical analysis on the galaxy sample as a 
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Table 3. Substructure analysis. Columns 6 is the highest value of columns 3,4,5 and gives the final likelihood for the presence 
of substructure 
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whole could be started. We chose the most straightforward 
methods, using absolute magnitudes M and velocities nor- 
malized by cluster mean velocity and velocity dispersion, 
i.e. : 



Vi 



(v) 



(1) 



SO-galaxies cannot be excluded by an F-test (17 % like- 
lihood). The same is true for the difference between SO 
and S-galaxies (23 % likelihood of same underlying distri- 
bution), while the difference between E and S-galaxies is 
significant at a level of 3 %. 



This implies that Vi has a mean value of and a stan- 
dard deviation of 1 in each cluster. Normalized velocity 
dispersions 



N(N — 1) 



(2) 



will be used in the following as kincmatical indicators 
for different samples. Only galaxies brighter than M = - 
igmag nave b een included in the following analysis, which 
roughly corresponds to the completeness limit of the 
present dataset. Choosing only objects brighter than 
M = -19 also excludes dwarf galaxies, which could con- 
taminate fainter samples (Binggeli etal. 1988). 

4-1. Dependence on galaxy types 

First, we looked at differences in the kinematical be- 
haviour of galaxies depending on their morphology. It can 
be seen in Figure 1 that there is a continuous trend of the 
velocity dispersion to increase from early to late galaxy 
types. Velocity dispersions have been computed with a bi- 
weight estimator of scale (Beers etal. 1990), which has 
shown to be superior when only few objects are involved. 
Error bars come from a bootstrapping calculation with 
1000 iterations. Between E and S galaxy types a rise in the 
velocity dispersion of 30 % can be observed. The hypoth- 
esis that E-galaxies have the same velocity dispersion as 




Fig. 1. Velocity dispersion for galaxies of different morpholog- 
ical classes. 
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4-2. Dependence on galaxy luminosity 

In the present sample of galaxy clusters there are clear 
signs of velocity dispersion dependence on absolute mag- 
nitude for galaxies brighter than M w — 22 maff , as can 
be seen in Fig. 2. Again, velocity dispersions have been 
computed with a biweight estimator of scale and error 
bars come from bootstrapping with 1000 repetitions. Bin 
limits were set every 0.5 mag between M = — 23.5 mas and 
M = — 19.0 maff , taking the biweight mean value < M > in 
each bin as the x-position instead of the bin center. Fur- 
ther dividing the sample into galaxies of different types 
reveals that mainly the most luminous early type galaxies 
are responsible for the lower velocity dispersion. However, 
there is a general tendency of ellipticals to have lower dis- 
persions than spirals also at the faint end. Former effect 
cannot be due exclusively to the existence of D /cD galax- 
ies residing in the bottom of the potential well, because 
only 4 of the clusters are given a Bautz-Morgan type I 
(Abcll etal. 1989). Another hint is that there are no E- 
galaxies with normalized velocities larger than 1.2 down to 
— 20.6 maff , while 6 or 7 would be expected from a normal 
distribution with a w = 1. No clear indication about the 
kinematical status of SO galaxies was found. It is widely 
known that there is considerable danger of confusion while 
classifying SO galaxies as an intermediate class between E 
and S (Bender 1992). This uncertainty in the morpho- 
logical classification, together with the low number of SO 
galaxies brighter than M = —22 mag in our sample, makes 
it difficult to state about the presence of luminosity segre- 
gation in this class. It seems that the brightest SO galaxies 
lie very close to the kinematical center of the cluster, indi- 
cating that these galaxies have been residing in the centers 
of clusters for long periods of time. On the other side, no 
signs of velocity dispersion changes with luminosity can 
be seen for E or SO galaxies fainter than M = — 21.5 mofl . 
This can be explained by the fact that the time-scale for 
dynamical friction exceeds the Hubble time when galaxies 
of luminosity lower than L* are involved (Sarazin 1988). 
As can be seen in Fig. 2, there is a deviating galaxy with 
M w -23.2 ma f and a w w 1.6, which was classified as S 
by Nilson (1973). Inspection of the corresponding photo- 
graphic plate reveals that the object is of peculiar nature, 
possibly interacting with its neighbours. Its morphologi- 
cal type is given as uncertain in several other catalogues, 
ranging between SO and Scd. 



5. Discussion 

We analysed the kinematics of the core regions of 12 
nearby galaxy clusters. The homogeneous redshift sample 
was nearly complete down to faint limits (< 18-19 ma9 ) 
and was supplemented by magnitudes and rough morpho- 
logical types for most of the galaxies. An analysis of sub- 
structure using velocity data revealed that 50 % of the 
cluster cores harbour significant substructure, thus con- 
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Fig. 2. Velocity dispersions for galaxies of different absolute 
magnitudes. 



firming that many clusters are not even relaxed in their 
inner regions, where effects like mass segregation and in- 
falling groups of galaxies might be disturbing the virial- 
ization process. 

After having merged the data to a sample of normalized 
galaxy velocities, magnitudes and types, we looked for 
type and type/luminosity segregation in velocity space. 
Previous findings (Binggeli etal. 1987; Sodre etal. 1989) 
about early type galaxies having lower velocity disper- 
sions than late types are confirmed by the present anal- 
ysis. Zabludoff & Franx (1993) found no such relation, 
on the opposite they claimed deviations in the velocity 
means between different types, concluding that there must 
be groups of spirals falling onto the cluster main body 
and distorting the distribution of velocities. Their findings 
should be considered complementary to ours, because of 
the different scale observed (R < 0.5 h _1 Mpc versus R < 
1.5 r 1 Mpc). 

Luminosity segregation in velocity space is also present, 
qualitatively and quantitatively in agreement with the 
findings of Biviano etal. (1992), who were using a larger, 
but more heterogeneous data sample. Moreover, there is 
a link between type and luminosity segregation. Only the 
brightest E and, possibly, SO-galaxies (M < -22 mag ) show 
clear signs of the phenomenon of luminosity segregation in 
velocity, which is probably related to two-body relaxation 
effects. These galaxies are responsible for the differences 
in velocity dispersion between early and late galaxy types, 
plausibly representing the fraction of galaxy population in 
clusters which have undergone significant late dynamical 
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evolution. On the other side, S-galaxies show no sign of 
luminosity segregation, as is expected from objects that 
are still infalling onto the cluster main body and are pre- 
sumably crossing the core region for the first time. 
It remains uncertain whether the effect is due to dynam- 
ical friction or comes from the fact that galaxies which 
are the result of repeated mergers will tend to have ve- 
locities closer to the cluster mean. In both cases, strong 
evolution in the dynamics of the oldest cluster population 
can be seen. Capelato etal. (1981) found several hints for 
the presence of partial equipartition of galaxy kinetic en- 
ergies in clusters. If the M /L ratio is close to constant for 
E-galaxies, then our findings support their view that only 
the most massive galaxies had had enough time to slow 
down. 
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